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Graphene is at the center of a significant research effort. Near-ballistic transport at room tem- 
perature and high mobility make it a potential material for nanoelectronics. Its electronic and 
mechanical properties are also ideal for micro and nanomechanical systems, thin-film transistors 
and transparent and conductive composites and electrodes. Here we exploit the optoelectronic 
properties of graphene to realize an ultrafast laser. A graphene-polymer composite is fabricated 
using wet-chemistry techniques. Pauli blocking following intense illumination results in saturable 
absorption, independent of wavelength. This is used to passively mode-lock an Erbium-doped fibre 
laser working at 1559nm, with a 5.24 nm spectral bandwidth and ~460 fs pulse duration, paving 
the way to graphene-based photonics. 



Ultrafast laser sources have many potential ap- 
plications, ranging from basic research and metrol- 
ogy to telecommunications, medicine and materials 
processingi*^!^. The majority of ultrashort lasers em- 
ploy a mode-locking technique, whereby a non-linear op- 
tical element - called saturable absorber - turns the laser 
continuous wave output into a train of ultrashort opti- 
cal pulses. Semiconductor Saturable Absorber Mirrors 
(SESAMs) currently dominate passive mode-locking^iii^. 
However, these have a narrow tuning range (tens of 
nm), and require complex fabrication and packaging^. 
A simpler and cost-effective alternative relies on Sin- 
gle Wall Carbon Nanotubes (SWNTs)^iii^, where the 
working wavelength is defined by choosing the SWNT 
diameter {i.e. bandgap)^. Tunability is possible by 
combining SWNTs with a wide diameter distribution^. 
However, when operating at a particular wavelength, 
the SWNTs not in resonance are not used, and con- 
tribute unwanted insertion losses, compromising device- 
performance. Novel nonlinear materials with broadband 
absorption are therefore required for wideband, tunable 
operation. 

The linear dispersion of the Dirac electrons in graphene 
offers the ideal solutionS.: for any excitation energy there 
will always be an electron-hole pair in resonance. Due to 
the ultrafast carrier dynamicsi^iiiii^ and large absorption 
of incident light per layer (ai = 2.3% — i^), graphene 
should behave as a fast saturable absorber over a wide 
spectral range. Compared to SESAMs and SWNTs, 
graphene saturable absorbers would not need bandgap 
engineering or chirality / diameter control to optimize de- 
vice performance. Here, we demonstrate an ultrafast fi- 
bre laser mode-locked at ^ 1.5^m, the most common op- 
tical telecommunications wavelength, using single layer 
graphene (SLG) and few layer graphene (FLG) fiakes. 

For ease of integration and stability, we incorporate 
the fiakes into a host polymer matrix. For the compos- 
ite preparation, we process graphite in water without 
functionalizationi^. This allows us to retain the elec- 
tronic structure of pristine graphene in the resulting ex- 
foHated SLG and FLG fiakes^^. We employ bile salts in 
order to obtain a stable, higher concentration of SLG 



and FLG than in previous non-aqueous dispersions^^. 
These amphiphilic molecules, with a hydrophobic and a 
hydrophilic side^^, disperse graphene in water by physical 
adsorption on its surface. In contrast to linear chain sur- 
factants, e.g. sodium dodecylbenzene sulfonate (SDBS) 
widely used for SWNTs, the fiat molecular structure of 
sodium deoxycholate (SDC), used here, efficiently dis- 
perses SLG, FLG and graphitic fiakes. Polyvinyl alco- 
hol (PVA) is chosen as host polymer for its mechani- 
cal properties and solvent compatibility. The dispersions 
are mixed with PVA via ultrasonication. The solvent 
is then evaporated to obtain free-standing composites 
~50^m thick. This simple, wet-chemistry approach is 
scalable and, more importantly, allows easy integration 
into a range of photonic systems (see Methods). 

We first characterize the dispersions obtained from 
ultrasonication, followed by centrifugation to remove 
large residual graphitic particles. Since we need to in- 
tegrate our composites in transmissive mode in a fi- 
bre laser cavity, fiakes with dimensions > l//m are re- 
moved to avoid optical scattering losses^. FigUJa) shows 
a photograph of a dispersion after centrifugation and 
its absorption spectrum. This is mostly featureless, as 
expectedi^>i^»i^>i^ii^. The peak in the UV region is a sig- 
nature of the TT plasmon2£. Using the experimentally de- 
rived absorption coefficient of 1390 L.g~s.m~s at 660nm 
(1.88 eV)^, we estimate ~0.08g/L of graphitic material 
in the centrifuged dispersion. 

The decanted dispersion mostly contains sub- 
micrometer fiakes. A typical Raman spectrum of such 
fiakes measured at 514. 5nm is plotted in Fig[T][b). 
Besides the G and 2D peaks, this has significant D and 
D' intensities, and the combination mode D+D'^2950 
cm~s. The G peak corresponds to the Fi2g phonon at 
the Brillouin zone centre. The D peak is due to the 
breathing modes of sp^ rings and requires a defect for 
its activation by double resonance (DR)2ii22ii2^. The 
2D peak is the second order of the D peak. This is a 
single band in monolayer graphene, whereas it splits 
in four in bi-layer graphene, refiecting the evolution 
of the band structure^i. The 2D peak is always seen, 
even when no D peak is present, since no defects are 
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Figure 1: (a) Absorption spectrum of a stable, 10% diluted, graphene dispersion. The insert shows a photograph of the 
undiluted dispersion. (b)Raman spectrum measured at 514. 5nm for a representative flake. (c) Absorption spectra of graphene- 
PVA composite and a reference PVA film. Insert: micrograph of a graphene-PVA composite. (d)Typical transmission of the 
composite as a function of average pump power for six excitation wavelengths. 



required for the activation of two phonons with the same 
momentum, one backscattering from the other. DR can 
also happen as intra-valley process, i.e. connecting two 
points belonging to the same cone around K or K'. This 
gives rise to the D' peak. The 2D' is the second order of 
the D' peak. The very large intensity of the D peak in 
Fig. lb is not due to the presence of a large amount of 
structural defects, otherwise it would be much broader, 
and G, D' would merge in a single band^^. We rather 
assign it to the edges of the sub-micrometer flakes we 
produce. We note that the 2D band, although broader 
than in pristine graphene^i, is still fitted by a single 
lorentzian lineshape. This implies that, even if the 
flakes are multi-layers, they are electronically almost 
decoupled and behave, to a flrst approximation, like a 
collection of single layers, retaining the Dirac fermions 
linear dispersion^. 

The micrograph of the composite shown in the insert 
of Fig[T][c) confirms the homogeneous distribution of the 
flakes. This is vital to reduce Mie scattering, that could 



be caused by flakes of dimensions comparable to 1.5/im, 
the device operation wavelengtbi^. The absorption spec- 
tra of the graphene-PVA composite and a reference PVA 
fllm are presented in Figdjc). The absorption of the 
flakes is featureless with the characteristic UV plasmon 
peak, while the host polymer only contributes a small 
background for longer wavelengths. 

FiglHd) plots the measured transmission as a function 
of average pump power at six different wavelengths (using 
a probe laser with pulse width 580 fs, as detailed in Meth- 
ods). At a relatively low input power level, the transmis- 
sion, r, is almost independent of pump power. However, 
T increases by At = 1.3% due to absorption saturation 
when the incident average power is raised to 5.35 mW 
(corresponding to a peak power density of 266 MW/cm^) 
at 1558 nm. Further increase in transmission is feasible, 
but limited to our maximum available pump intensity. 
The saturable absorption of the composite is clear for 
all pump wavelengths. This indicates that our composite 
can be used over a broad spectral range, unlike SESAMs^. 
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Figure 2: Schematic illustration of photoexcited electron ki- 
netics in graphen. The intensity of the color shading rep- 
resents the electron population, in the case of (a) efficient 
interband relaxation and (b) inefficient interband relaxation. 



The transmission change is similar to SESAMs^. The 
non-saturable insertion loss (34.3%) is larger than that 
of SESAMs^, but comparable to SWNTs^. Note that, 
for fibre lasers with relatively large single round-trip gain 
coefficient, such non-saturable losses are tolerable^. Fur- 
ther decrease in non-saturable insertion loss is expected 
when the device is completely saturated. To estimate the 
non-saturable loss due to coupling between the two fibre 
connectors, a reference ~ bQ^m PVA composite is used 
in place of the graphene-PVA composite. The transmis- 
sion of the packaged pure PVA is tq = 82.6%. Given 
the pulse repetition rate of 38.83 MHz, we can estimate 
the density of photons absorbed per pulse at our maxi- 
mum pump power to be 2.3 x 10^'* cm~^. The average 
number of layers N participating in the absorption is es- 
timated from r = ro(l — Qfi)^, where ai = 2.3% is the 
absorption per layer, yielding A*" ~ 11 (see Methods). 
Thus, the photon density absorbed per pulse per layer is 
2.2 X lO^^cm"^, and Ar = 1.3% translates into a relative 
change in the absorption per layer Aai/ai = —8.2%. 

To understand the origin of the observed saturation 
and give its theoretical estimate, let us discuss the pho- 
toexcited carrier kinetics in graphene (Fig. [2]) . Optical 
interband excitation in SLG and FLO by an ultrashort 
optical pulse produces a non-equilibrium carrier popu- 
lation in the valence and conduction bands. In time- 
resolved experiment a^°i^^ two relaxation time scales are 
typically seen. A fast one, ^ 100 fs, usually associated 
with carrier-carrier intraband collisions and phonon emis- 
sion, and a slow one, on a picosecond scale, correspond- 
ing to electron interband relaxation and cooling of hot 
phonon a^^i^^ . Several points should be noted here, (i) In 
our experiment, SLG and FLG flakes are incorporated in 
a polymer matrix, so we assume graphene phonons to ef- 
ficiently give away their energy to the matrix and neglect 
hot phonon effects, (ii) The rate of carrier-carrier colli- 
sions grows with carrier density, i. e., pump power. We 
assume the coUisions to be very fast.(iii) Carrier-carrier 
collisions preserve the total electronic energy and thus 



cannot change the effective electronic temperature. 

A quantitative treatment of the electron-hole dynam- 
ics would require the solution of the kinetic equation for 
electron and hole distribution functions /e (p) and fh (p) , 
p being the momentum counted from the Dirac point. 
Here we present an estimate. Let us assume that relax- 
ation times are shorter than the pulse duration, so during 
the pulse the electrons reach a stationary state, that colli- 
sions put electrons and holes in thermal equilibrium at an 
effective temperature T^ff, and take undoped samples: 



/e(p) = A(P) 



1 



(1) 



where v is the velocity of the Dirac electrons in SLG. 
The populations determine the electron and hole densi- 
ties ne,h and the total energy density £ (counted from 
the energy of the un-doped sample at zero temperature): 



ne,h = 4 



£ 



d2p 



/e,/i(p), 



(2) 



^'p[/c(p)+/h(p)] 



the factor 4 accounts for valley and spin degeneracy. 
The populations also determine the reduction of photon 
absorption per graphene layer for a given laser energy 
El, due to Pauli blocking, by a factor 1 + Aai/ai = 
[1 — /e(p)][l — //i(p)]. In our measurements, vp = El/ 2 = 
0.4 eV and 1 + Aai/ai ~ 0.92, implying Tg// 0.13 eV. 

To obtain a theoretical estimate for T^f / , let us assume 
that during the pulse electrons and holes are injected at 
the energy = £'l/2 = 0.4 eV at a constant rate. 
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This corresponds to pumping the energy at a rate: 
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2 J,„e,„ = 3.0 • 10^^ cm-^ps-^eV. (4) 
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The electronic energy density £ can be decreased by 
phonon emission. Neglecting hot phonon effects, for the 
Fermi-Dirac distribution Jl]) the cooHng rate can be cal- 
culated as (see Methods): 
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(5) 
(6) 



The first two terms in Eq. ^ correspond to emission 
of optical phonons with wave vectors near F and K, re- 
spectively. Their frequencies lOrTi^K are assumed to be 
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independent of wave vector (Einstein model). The di- 
mensionless coupling constants Ar and \k are defined 
in Ref. [IS- The third term in Eq. ([5]) corresponds to 
emission of longitudinal acoustic phonons, with frequency 
taken proportional to their momentum (Debye model). 
The length Lac in Eq. ^ is expressed in terms of de- 
formation potential, Do, carbon atom mass, Af, and 
unit cell area, ^u.c, as Lac = D^A^^cJi'^Mv'^). For 
Dq = 10 eV, ^u.c. = 5.24 A^ M = 2.00 x lO^^s g ^ 
2.88 X 10^ (eV-A^)-\ v = 10*cm/s = 7eV-A, we obtain 
Lac = 0.26 X 10-3 A. We take cur = 0.20 eV, Ar = 0.03, 
luk = 0.15 eV, Xk — 0.1—. These give the time of opti- 
cal phonon emission by the photoexcited electron to be 
[(Ar/2)(L;L/2 - a;r) + {Xk/2){El/2 - cokT' = 40fs21. 
From the thermal balance equation. 



Pump laser 
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we obtain a theoretical estimate Tg// = 0.20eV, the dom- 
inant contribution to cooling coming from optical phonon 
emission. Given this value of Te//, and using the Fermi- 
Dirac distribution ((l]), we obtain that the absorption of 
photons with energy ~ 0.8 eV should be reduced by 
up to 1 Aai/ai = (1 - /e)(l - fh) = 0.78. This cor- 
responds to Aq!i/q;i = —22%, while the measured one 
is Aai/ai = —8.2%. Thus, having assumed efficient 
carrier-carrier relaxation (both intraband and interband) 
and efficient phonon cooling, with the main bottleneck 
being the transfer of energy from the electrons to the 
graphene phonons, we obtain a higher saturation level 
than the measured one. On one hand, this points to 
other mechanisms contributing to electron cooling (e. g., 
interaction with the polymer matrix). On the other hand, 
more importantly, it suggests that sample optimization 
could lead to enhanced saturable absorption and, thus, 
much better performance and lower power consumption. 

The above estimates are made under the assumption of 
efficient relaxation due to carrier-carrier collisions, both 
intraband and interband. However, in the case of lin- 
ear dispersions near the Dirac point (i.e., for SLG), pair 
carrier-carrier coUisions cannot lead to interband relax- 
ation, thereby conserving the total number of electrons 
and holes separately. A three-particle collision is required 
to move an electron from the conduction to the valence 
ban d^^i^" . Interband relaxation by phonon emission can 
occur only if the electron and hole energies are close to 
the Dirac point (within the phonon energy). Note that 
for graphite and graphite flakes the situation is com- 
pletely different: the dispersion near the Dirac point is 
quadratic, and pair carrier collisions can lead to inter- 
band relaxation. Thus, in principle, decoupled graphene 
monolayers (as in our present experiment, as shown by 
the FLG Raman spectra) will provide the highest sat- 
urable absorption for a given amount of material. To pre- 
cisely estimate such effect, three-particle collisions would 
have to be considered. In our case, interband relax- 
ation due to phonon emission seems insufficient to bal- 
ance the pumping of electrons and holes at the rate given 
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Figure 3: Graphene mode-locked fibre laser and mode-locker 
assembly. ISO: isolator; WDM: Wavelength division multi- 
plexer; PC: Polarization controller;EDF: Erbium doped fibre. 



by Eq. ^ (see Methods). Thus, three-particle collisions 
do play a role. This is not surprising, since the dimen- 
sionless Coulomb coupling constant Ts is not very small, 
and the effective temperature, Tg// = 0.13 eV, is quite 
high for the corresponding time scale, (rfTe//)"^, to be 
in the femtosecond range. 

We now use our composite to build and test an ul- 
trafast laser working at the main telecommunications 
window. The mode-locker is assembled by sandwich- 
ing the graphene-PVA between two fibre connectors with 
a fibre adapter, as schematized in FigO A 0.8m Er- 
bium doped fibre (EDF) is used as the gain medium. 
It is pumped by a 980nm diode laser via a wavelength- 
division-multiplexer (WDM). An isolator (ISO) is placed 
after the gain fibre to maintain unidirectional operation. 
A polarization controller (PC) optimizes mode-locking. 
A 20/80 coupler is used. The total cavity length is 
10.54 m. The input mode diameter is ^^10 yum. 

The threshold pump power for continuous wave las- 
ing is ~10mW. When this is increased to ~27mW, sta- 
ble mode-locking can be initiated by introducing a dis- 
turbance to the intracavity fibre. Once stable output 
is achieved, no further polarization controller adjust- 
ment is needed. It is always possible to decrease the 
pump power to ~22mW while maintaining mode-locking. 
When mode-locked, the laser produces a pulse train at a 
rate of 19.9MHz. 

FiglUa) shows a typical output spectrum, with cen- 
tral wavelength ~1559nm. The full width at half max- 
imum (FWHM) bandwidth is 5.24nm. The sidebands 
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Figure 4: Mode-locked pulses characteristics. (a) Output spectrum; Spectral resolution O.lnm.(b) Autocorrelation trace of output 
pulses;Delay resolution~20fs.(c)Oscilloscope trace. (d) Wideband RF spectrum up to 0.8GHz. (e)RF spectrum measured around 
the fundamental repetition rate, /i = 19.9 MHz.(f)RF spectrum measured around the tenth harmonic of the repetition rate, 
/io=199MHz. In (e,f), the measurement span is 2KHz with 1-Hz resolution bandwidth. The red trace depicts the background 
when the laser is switched off. 



at 1548.6 and 1570nm, are typical of soliton-like pulse 
formation^, resulting from intracavity periodical per- 
turbations^. FiglUb) is a second harmonic generation 
(SHG) autocorrelation trace of the mode-locked pulses 
after a 50cm single mode fibre lead from the coupler, with 
FWHM=713 fs. Assuming a sech^ temporal profile, the 
deconvolution gives a pulse duration of 464fs. The time- 
bandwidth product (TBP) is 0.3. Due to the autocorrela- 
tor delay resolution (20fs,5%), there is a minor deviation 
from 0.315, expected for transform-limited sech^ pulses, 
having the shortest duration for a given spectral width. 



This confirms soliton-like operation^^. Figure |4][c) is the 
output pulse train, with a period 50.3ns, as expected 
from the fibre cavity length. To study the operation sta- 
bility, we measure the radio-frequency (RF) spectrum. 
We first perform wide span frequency spectrum mea- 
surements up to 800MHz, as in Fig|4ljd). This shows no 
significant spectral modulation, implying no Q-switching 
instabilities^. Figsl4Ue,f) are the RF spectrum around 
the fundamental repetition rate and its tenth harmonic. 
A >70dB peak-to-background ratio (10^ contrast) is ob- 
served, indicating good mode-locking stability^. 
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Our graphene-based ultrafast laser, harnessing the 
wideband optical nonlinearity of graphene, with no need 
of bandgap engineering, extends the practical potential 
of this novel material from nanoelectronics to optoelec- 
tronics and integrated photonics. 

We acknowledge funding from EPSRC GR/S97613/01, 
EP/E500935/1, the European Research Council grant 
NANOPOTS, a Royal Society Brian Mercer Award for 
Innovation, the Isaac Newton Trust and the Cambridge 
Integrated Knowledge Centre in Advanced Manufactur- 
ing Technology for Photonics and Electronics. 



Methods 

A. Sample Preparation and Characterization 

Dispersions are prepared by ultrasonicating (Bran- 
son 450A, 180W power, 8-10°C) 120mg graphite (Sigma 
Aldrich) for 2hrs in 10ml deionized (DI) water, with 90mg 
Sodium deoxycholate (SDC) bile salt. The resulting dis- 
persion is then centrifuged at 5 krpm (2500^) using a TH- 
641 swing rotor in a Sorvall WX-100 Ultra centrifuge for 
Ihr at 18°C. The PVA-graphene composite is prepared by 
mixing 4ml of the centrifuged dispersion with 2ml 15wt% 
aqueous PVA (Wako chemicals) solution in a Hauschild 
Speedmixer (DAC150 FVZ-K) for 5min at 3krpm. The 
mixture is then vacuum evaporated at 20°C, resulting in 
a ~50/im thick composite. 

For UV-Vis-NIR absorption spectroscopy, 1 part of 
centrifuged dispersion with 10 parts of 0.9wt% SDC- 
water solution are mixed. The absorption is background 
subtracted to account for solvent and surfactant. A 
Perkin-Elmer Lambda 950 UV-vis spectrometer is used. 
The centrifuged dispersions are diluted as before and 
drop-cast onto a Si wafer with 300nm thermally grown 
Si02. The drops are dried at 20°C. The samples are then 
rinsed in acetone for 5mins and water for an additional 
5 mins. This removes the surfactant. These samples are 
then used for Raman and SEM characterization. Raman 



J 



spectra are collected with a Renishaw 1000 spectrometer 
at 514. 5nm using a lOOx objective, and ~lmW power. 

For power-dependent transmission measurements, the 
graphene mode-locker assembly is coupled to a ~580fs 
source, tuneable from 1548 to 1568 nm. The typical spec- 
tral bandwidth at 1558nm is 2.8nm. This is achieved 
by filtering a homemade SWNT fibre laser (^400 fs, 
38.83MHz) using a 3-nm band-pass filter. The laser beam 
is then amplified by an Erbium-doped fibre amplifier and 
a 20% tap monitors the input power. Another 80% is 
used to pump the graphene mode-locker assembly. Two 
calibrated power-heads are programmed to read the in- 
put/output power simultaneously. The estimated effec- 
tive mode area on the composite, A^ff, is 78.54 ^m^, 
deduced from the ^ 10/im mode diameter. 

A spectrum analyzer (Anritsu MS9710B) and a SHG 
autocorrelator (APE Pulsecheck 50) are used to measure 
the output spectrum and pulse width. The RF spectrum 
is measured using a home-made lightwave converter con- 
nected to a spectrum analyzer (Anritsu MS2719B). 

For a pump wavelength of 1558 nm, the total number 
of pump photons per cm^ per pulse (Ttot) at 5.35 mW is 



5.35 mW 



38.83 MHz X Elx A^ff 



1.38 X 10 



i5cm-2. 



As discussed in the main text, the measured transmission 
of a packaged reference pure PVA device is tq — 82.6%. 
Considering the absorption of N graphene layers, 1 — (1 — 
Q!i)^, with ai — 2.3%, and Hnear coupHng losses, 1 — tq, 
to cause the non-saturated loss 1 — r (t = 65.7%), the 
total number of photons absorbed by the SLG and FLG 
fiakes per cm^ per pulse for our maximum pump in the 
power-dependent absorption measurements is (82.6% — 
65.69%) X Ttot = 2.3 x lO^'^cm'^. This can be estimated 
to be on average equivalent to the effect of ^11 graphene 
layers: 65.7/82.6 ~ (1 - 2.3%)i°-^. Then, the average 
number of photons absorbed per layer per cm^ per pulse, 
is - (2.3 X 10i'*cm-2)/li = 2.2 



B. Kinetic Equations 



The kinetic equations for the evolution of electron and hole distribution functions /e(p),//i(p) due to phonon 
emission can be written analogously to Ref. 27. Here the task is even simpler as (i) any spatial dependence is absent, 
(ii) due to frequent coUisions each photoexcited electron-hole pair is assumed to quickly forget its initial state, so the 
joint electron-hole distribution function /(p,p') factorizes into a product /e(p) //i(p')- For fe,h{p) we have: 



dfejp) 

dt 



(8) 



dfkip) _ f v'd'v 



dt 



J 7^ ^^^P'' P^ ^'^^P'^t^ ^'^^P^l " ^^P' P') -^"^P^I^ " ■^'^^P')]^ J 7^ ^^P' P'^ -^^^P) ^''^P'^- 



(9) 

The phonons are assumed to be at zero temperature. The interband relaxation kernel i?(p, p') is a sum of three terms, 
corresponding to emission of (i) optical phonons with wave vectors near K, (ii) optical phonons with wave vectors 
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near F, and (iii) longitudinal acoustic phonons: 

i?(p, p') = Rk{p, p') + Rr{p, p') + i?ac(p, p'), 
Rk{p, p') - 2Xk sin^ ^fR-JfEL T,s{vp - vp' - ujk), 



Rt{p,p') = XrTrS{vp~ vp' - wr), 



-Rac(P,p') 



2 (pp - V?p' 



TrS{vp — vp' — Vs\p — p' 



The interband kernel contains only the optical phonon contribution: 

i?(p,p') = i?K(p,p')+i?r(p,p'), 

Rk{p, p') = 2Xk cos^ ~ "^P' 7r(5(z;p + V - c^i^), 

^r(p, p') = Ar 7r(5(wp + wp' - tor)- 



(10) 

(11) 
(12) 
(13) 

(14) 
(15) 
(16) 



Here we use the Einstein model for the optical phonon dispersion (momentum- independent frequencies u!K,uJr) and 
the Debye model for the acoustic phonon dispersion (the frequency of the phonon with wave vector q is Ws|q|, where 
Vs <^ V is the sound velocity). We introduce the dimensionless electron-phonon coupling constants Xk^^f for the 
optical phonons, while the analogous coupling constant for the acoustic phonons has the dimensionality of length, 
since the coupling is proportional to the phonon wave vector: 



Mv^ujr 



Xk = 



27a' 



^0 



27a' 



Mv^v^ 



(17) 



Here Fk and Fr are the force constants measuring the change in the electron energy with lattice displacement along 
the corresponding normal mode (in the nearest-neighbor tight-binding model Fr = Fk = 3{dto/da), where dto/da is 
the derivative of the nearest-neighbor electronic coupling matrix element to with respect to the bond length a) , Dq is 
the deformation potential, M is the mass of the carbon atom, V27a^/4 is the area per carbon atom (half of the unit 
cell area). These definitions coincide with those of Ref. 27. 

For a rigorous treatment of the problem, this kinetic equation should be supplemented with a source term describing 
the optical pumping of electron and hole populations, as well as electron-electron collision terms for both pair and 
triple collisions, and then solved. For an estimate, we assume the electron-electron relaxation to be efficient enough 
to quickly establish the Fermi-Dirac distribution (U) and calculate the cooling rate due to the phonon emission. 
Combining the definition of £ [Eq. (2)] with the kinetic equation, we get: 



ph 



d^p fv^d^p' 



{/e(p) [1 - /e(p') + Mp) [1 - h{p% X 



(27r)2 J (27r) 

X [ujKRKip,p') +i^rRrip,p') + Vs\p - p'| i?ac(p,p')] + 

d^p f v'^d'^p' \ ~ ~ 

T^r^ / /n \2 fe{p)fh{p') w/fi?i^(p,p') +^ri?r(p,p') 
(27r)^ J (27r)^ L 



(18) 



Substituting here the Fermi-Dirac distribution ((T|) and calculating the integrals in the approximation Vs 'ti v, we 
arrive at Eq. (5), with Lac = ^ac^s/w. 

The rate of interband relaxation due to phonon emission is calculated analogously. 



dup 



dt 



ph 



d'p f v'^d'^p' 



(27r)2 J (27r) 



/e(p)A(p') i?K(p,p') + ^r(p,p') 



< 



(19) 



The upper limit on the relaxation rate is obtained by setting /e = fh = 1- We get 9 x 10^^ cm ^ps ^, which is four 
times smaller than the pumping rate in Eq. (3), thus valida ting our assumptions. 
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